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Abstract 

Using a quantum mechanical model, the exact energy eigenstates for two-particle 
CN ' two-channel scattering are studied in a cubic box with periodic boundary conditions 

in all three directions. A relation between the exact energy eigenvalue in the box and 
the two-channel S'-matrix elements in the continuum is obtained. This result can be 
viewed as a generalization of the well-known Liischer's formula which establishes a 
■ similar relation in elastic scattering. 
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^ ■ 1 Introduction 



Scattering experiments play an important role in the study of interactions 
among particles. In these experiments, scattering cross sections are measured. 
By a partial wave analysis, one obtains the experimental results on particle- 
particle scattering in terms scattering phase shifts in channels of definite 
quantum numbers. In the case of strong interaction, experimental results on 
hadron-hadron scattering phase shifts are available in the literature [1,2,3,4,5]. 
On the theoretical side, Quantum Chromodynamics (QCD) is believed to 
be the underlying theory of strong interactions. However, due to its non- 
perturbative nature, low-energy hadron-hadron scattering should be stud- 
ied with a non-perturbative method. Lattice QCD provides a genuine non- 
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perturbative method which can tackle these problems in principle, using nu- 
merical simulations. In a typical lattice calculation, energy eigenvalues of two- 
particle states with definite symmetry can be obtained by measuring appro- 
priate correlation functions. Therefore, it would be desirable to relate these 
energy eigenvalues which are available through lattice calculations to the scat- 
tering phases which are obtained in the scattering experiment. This was ac- 
complished in a series of papers by Luscher [6,7,8,9] for a cubic box topology. 
In these references, especially Ref. [8], Luscher found a non-perturbative re- 
lation of the energy of a two-particle state in a cubic box (a torus) with the 
corresponding elastic scattering phases of the two particles in the continuum. 
This formula, now known as Liischer's formula, has been utilized in a number 
of applications, e.g. linear sigma model in the broken phase [10], and also in 
quenched QCD [11,12,13,14,15,16,17,18]. Due to limited numerical computa- 
tional power, the s-wave scattering length, which is related to the scattering 
phase shift at vanishing relative three momentum, is mostly studied in hadron 
scattering using quenched approximation. CP-PACS collaboration calculated 
the scattering phases at non-zero momenta in pion-pion s-wave scattering in 
the 1 = 2 channel [17] using quenched Wilson fermions and recently also in 
two flavor full QCD [19]. 

For hadron scattering at low energies, usually the elastic scattering is dominant 
since the inelastic channels are not opened. However, when the energy of 
the scattering process exceeds some threshold, inelastic scattering starts to 
contribute and the scattering of the particles cannot be described by single 
channel elastic scattering anymore. In the case of pion-pion scattering, for 
example, the scattering process is elastic below the four pion and the two kaon 
threshold. If the center of mass energy exceeds the four pion threshold, inelastic 
effects starts to contribute. The inelastic effects become very important when 
the energy is getting close to the two kaon threshold. At this point, pion- 
pion can be scattered into kaon-kaon pair final state. Although the four pion 
threshold is in fact below the two kaon threshold, four pion final state will not 
contribute significantly due to its weak coupling to the two pion initial state. 
The fact that four pion states are coupled to the two pion state weakly in the 
low-energy limit is seen from the QCD chiral lagrangian. A six pion vertex 
in this lagrangian involves derivative coupling which is vanishing in the low- 
energy limit. Experimental investigations show that the contribution of four 
pion states would only make substantial contribution when the energy is well 
over lGeV. Therefore, if the energy is not much higher than lGeV, pion-pion 
scattering can be approximated rather well by a two-channel model. It is then 
interesting to study the relation between the multi-channel two particle states 
and the scattering phases, just as what we have done in the single channel 
case. 

In this paper, we establish a relation between the energy of a two particle state 
in a finite cubic box and the scattering matrix parameters. It is a generalization 
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of the famous Liischer formula to the multi-channel situation. This relation is 
non-perturbative in nature and it is derived in a quantum mechanical model 
of two-channel scattering. The result can also be generalized to the case of 
asymmetric box. Further generalization to the case of massive field theory is 
under consideration. 



2 The quantum mechanical model to two channel scattering 



In this paper, we study a quantum mechanical model of two-channel scatter- 
ing. Generalization to more channels can be done similarly. The model under 
investigation has the following Hamiltonian: 



E T-^ 2 ) \^{r)V 2 {r) J 



We will also use the notation: 



V(r) 



V x {r) A(r) 
A*(r) V 2 (r) 



(2) 



to denote the matrix valued potential whose matrix elements Vi(r), V 2 (r) and 
A(r) all vanish for r > R. E T > designates a positive threshold energy of 
the second channel. That is to say if the center of mass energy for the scatter 
process is less than E T , there will be no asymptotic scattering states in the 
second channel. 

Energy eigenstates of the Hamiltonian (1) with energy E can be decomposed 
into spherical harmonics: 

*(r) = £yu*)|**^| . P) 



l,m 



The radial wave-functions ipi-i m (r) with i — 1,2 satisfy the radial Schrodinger 
equation: 



+ - 3 +E-V l {r) 



2m\ \ dr 2 r dr r 2 
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d 2 2d 
2m 2 \ dr 2 r dr 




+ E' - V 2 {r) 



ip2(r) =A*(r)ipi(r) 



(4) 



where E' = E — E?- 

About the coupled differential equations (4), the following statement can be 
proven. 

Theorem 1 If the matrix valued potential (2) is such that every matrix ele- 
ment of r 2 V(r) is analytic around r = and that lim r ^ r 2 V(r) = 0, then the 
coupled differential equations (4) has two finite, linearly independent solutions 
near r = 0: u^j(r), with i = 1,2 designating different solutions and j = 1,2 
denoting different component of the solution. Moreover, these solutions can be 
chosen such that: 



u 



(5) 



Equations (4) has two further linear independent solutions vf.^r) that are 
unbounded near r = and satisfy: v\^(r) ~ r~ l ~ 1 5 i j. 

The proof of this theorem is quite similar to the proof in the single channel 
case. Details are provided in the appendix. 



3 Two channel scattering and the S matrix in infinite volume 



In this section, we briefly describe the quantum mechanical treatment of in- 
elastic scattering. We will concentrate on the two channel case, although the 
formalism can easily be generalized to more channels. 

At large r where the potential V(r) vanishes, the wave function of the scat- 
tering state can be chosen to have the following form: 



V 



e ikl ' r + /n(ki-f) 
f2i(k-r), 



i rn,2 e 



i k 2 T 



(6) 



This wave function has the property that in the remote past, it becomes an 
incident plane wave in the first channel with definite wave vector It is an 
eigenstate of the full Hamiltonian with energy: E = k^/(2m!). Similarly, if the 
energy E > -Er> 

one can also build an eigenstate of the Hamiltonian which in 
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the remote pase becomes an incident plane wave in the second channel: 
^ 2 )(r) r - 



' /i 2 (k 2 -f)J5^ 



Ve lk2 - r + / 22 (k 2 -r> 



/A ; 



(7) 



This state is also an eigenstate of the Hamiltonian with energy E = E T + 
k|/(2m 2 ). 

In partial wave analysis, one decompose the coefficients: appearing in the 
above construction into spherical harmonics: 



/n(ki-r) 
/ 2 i(ki-r) 



1 



2ih 
1 



£(2J + 1) sg(fci)-l fl(k!.f), 



^(2Z + l)^?(^)PKki-r), 



(8) 
(9) 



and similar expressions for / 12 and / 22 . With this we find the two scattering 
eigenstates in Eq. (6) and Eq. (7) can be expressed as: 



$W(r)^E^(^;(ki) 



lm 



1 



2ik\r 



^(Ogifcir _|_ /y+lg-ifcir 



\ 2i\/k\k2r 



°21 e 



*(2)( 



X 4 ^-(f)^(k 2 ) 

/m 



\ 2 jfc 2 T 



1 /mT^Wgifcir 
^(0 e ife 2 r _|_ ^_y+lg-ifc 2 r 



• (io) 



The two component wave-functions appearing in the above formulae are in 
fact the radial wave functions of the Schrodinger equation (4) in the large r 



region 



. 2 



(r) 



w\ >{r) 



i 



2ik-i r 



^(Ogifcir _|_ ^_y+l e -ikir 



1 



2i\/k\k2T y mi 



m 2 c(0„ifc 2 



1 / mi 

- ^i 2 e 



2 iV k\kir y m 2 



1 

\ 2 ife 2 r 



gWgifc 2 r _|_ ^ y+lg— ife 2 r 



11^ 



2 In fact, if we use the corresponding spherical Bessel's function ji and n;, we can 

3t an express 
when r — > oo. 



get an expression for to^(r) and iw^(r) for r > R. Eq. (11) is the asymptotic form 
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It is obvious that the two radial wave functions u>^(r) and wl 2 \r) are linearly 
independent. Therefore, according to the theorem stated in the previous sec- 
tion, they are linear superpositions of the general solutions: u^\r) and u[ 2 \r). 
The converse is also true. The radial wave functions u^\r) and u\ 2 \r) defined 
via Eq. (5) can be expressed as linear superpositions of the two radial wave 
functions in Eq. (11). In other words, there exists a non-singular 2x2 matrix 
C such that: 

™?\r) = E C^ir) , «»(r) = £ C^w^r) • (12) 
j j 



Another important physical property of the wave-function (10) is that the 
matrix elements which enter the expansion, namely form a 2 x 2 unitary 
matrix which is nothing but the S'-matrix in the subspace with orbital angular 
momentum /. This unitarity condition comes directly from the probability 
conservation law of quantum mechanics. We have, for example: 



|c(0,2 



'11 



|C(0|2 



'21 



1 



iq(0| 

l°12 I 



+ 



|C(0|2 



'22 



(13) 



In practice, if the theory has CP symmetry which is the case in QCD, the 
two-channel ^-matrix is usually parameterized as: 3 



S {l \E) 



I 



Vie 



2iS[ 



me 2iSi 



(14) 



where the real parameters: S[(E), S l 2 (E) and rji(E) are all functions of the 
energy E. We will assume in the following that the S'-matrix of the scattering 
problem has this form. 



4 Energy eigenfunctions on a torus 



Now we enclose the system discussed above in a cubic, periodic box with finite 
extension L in every spatial direction. The Schrodinger equation for the system 
now takes a similar form as in the infinite volume except that the potential 

3 For the potential models, we assume that the potential is invariant under time 
reversal and parity. Then, one has Sfi = Si*f*, where i* and /* denotes the time- 
reversed state of i and /. With this in mind, it is easily seen that eq. (14) is the 
most general 2x2 unitary matrix for the scattering matrix of spinless particles. 
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is periodically extended and the (two-component) eigenfunction has to satisfy 
the periodic boundary condition: 



[H + V L (r)]4>{r) = E^(r) , ^(r + Ln) = ^(r) . 



(15) 



where the free Hamiltonian is given by: 

/ ln2 n \ 



Ha 



2mi 







V 



E T -^j 



(16) 



and the periodically extended potential is: 
V L (r) = Y,V(r + Ln). 



(17) 



The eigenvalue equation (15) now has discrete spectrum and the corresponding 
eigenfunctions are smooth. 

It is convenient to partition the whole space into two regions. In the inner 
region, every point satisfies the condition: |r| < R, mod(L). In the outer 
region: 



Q = {r| : |r] > R ,mod(L)} . 



(18) 



Note that in the outer region Q, the interaction potential V^(r) = and the 
Schrodinger equation (15) reduces to two decoupled Helmholtz equations: 



(V 2 + A; 2 )^(r) = 0, i=l,2 , 



(19) 



where the energy eigenvalue E is given by: 

^1 



E 



2m\ 



7-1 ^2 

E T + - — 

2m 2 



(20) 



For energy < E < E?, k\ is real and hi takes purely imaginary values; for 
energy E < 0, both k± and /C2 are purely imaginary; for E > E?, which is the 
case of two channel scattering above the threshold, both k\ and k 2 are real 
numbers. 

It it easy to see that: 



*(r;£) = E 



Y lm (n) . 



(21) 



lm 



.i=l 
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solves the Schrodinger equation in the inner region for |r| < R. In the outer 
region Q, for a given value of energy E, the corresponding eigenfunction must 
be superposition of the free Schrodinger equation, which in the outer region 
decouples to two independent Helmholtz equations. Since there are two lin- 
ear independent radial wave functions, the eigenfunction must be some linear 
combination of the two: 



*(r;£) = £ 



lm . 



r) 



Y lm (n) , for r > R. , 



(22) 



with cf2 = J2i bimCij being non- vanishing coefficients. Note that when the 
system is enclosed in a finite periodic box, the exact energy eigenvalues become 
discrete. The degeneracy in the radial wave-function in general is then lifted. 
That is to say, for a given value of energy, there exists only one radial wave- 
function, unlike in the infinite volume where there are two such wave-functions 
for a given energy. 

On the other hand, in the region Q, the solution must be linear superposition 
of the singular periodic solutions of Helmholtz equation: 



tt(r;£) 



Elm v lmGlm( r ] kf) 
Elm Vlm G lm(r; fc£) 



(23) 



Combining Eq. (23) and Eq. (22), using the basic expansion of Gi m (r): 
(_)ijfci+i 



G lm (r;k 2 ) = 



47T 



Y lm (tt r )m(kr) + J2 Mim;l>m'Yy m >(tt T )ji>(kr) 
I'm' 



,(24) 



we arrive at the following set of linear equations: 



r (i) K w i n i r p),/^^o(o_ y H t,fc i +1 .,(i) M m 

c «mA D ll + L ) + c lrn y 12 _ ^ V l'm' Jvl l'm';lm > 

n: / x a , z'_i_i 

I (2) 



_ ?r (i) r c(0 

Jc lml°ll 

ij + Clm \ k.m, 21 47T V l'm'^l'm';lm 

I h n mn. n /» (_y/^ +1 (2) 
: 4tt ^ m ' 



r (2)/o(0 
C imA D 22 



_ ?V (2)/o(0 



, ( /. _ n _ <i) fe(o_(ii; (2) f25 x 

' lj ,Cta V*imi 21 " 4tt ^' (25) 

the symbol -M^/. Jm represents M.v m '-i m {kf) for i — 1,2, 
cplicit expression for for M.v m ';i m {kf ) are given in Ref. [8] 



In these equations, the symbol M$ r 
respectively. The explicit expression 
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which we quote here: 



M lm . js {k 2 ) = ]T 

I'm' 



(-)V-% m ,(l,g 2 

n 3/2 q l'+l 



7(2/ + l)(2/' + l)(2j + l) 



x 



I It A 




fir A 


v o 0) 




[mm'—sj 



(26) 



Here we have used the Wigner's 3j-symbols and q = kL/(2n). The zeta func- 
tion Zi m (s,q 2 ) is defined as: 



Zim{s,q 2 ) = J2 



yimjn) 
n 2 — q 2 ) s 



(27) 



According to this definition, the summation at the right-hand side of Eq. (27) 
is formally divergent for s = 1 and needs to be analytically continued. Follow- 
ing similar discussions as in Ref. [8], one could obtain a finite expression for 
the zeta function which is suitable for numerical evaluation [20,21]. From the 
analytically continued formula, it is obvious from the symmetry of O(Z) that, 
for I < 4, the only non-vanishing zeta functions at s = 1 are: Z 00 , and Z 40 . 



Eliminating the coefficients and from the set of equations (25) , one 
obtains a homogeneous linear equation for the coefficients and . In order 
to have non-trivial solutions for the coefficients an d c^, the corresponding 
matrix has to be singular. This condition then gives: 



(2) 



M {1)+ 

J 1 1' m'dm 



I'm' dm 



k 2 m 2 q{l') A ,f ( 2 ) 
o 21 JVl,,, 



k\m\ 



I'm' dm 



fcimi a /ft 1 ) 



k 2 m 2 



'12 



'M\. 



m';lm 



m'dm 



sg ] M\ 



(2)- 
m'dm 



(28) 



where the matrix elements M^y^'dm are defined as: 
r«± <i.i\ - n.2 



M^, lm (kf) = MZ'dm(kf) ± i^m'm , (29) 



with the parameter k 2 related to the exact energy eigenvalue via Eq. (20). 
Now if we further assume that the matrices: -i m {k 2 ) are non-singular, 

Eq. (28) may be expressed as: 

U l'm'dm ~ b ll d l'l d m'm \j ^^21 °l'l°m'm =Q 
\f^^ S ^l'l^m'm Vfm'dm ~ S 22^l'l^m'm 
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where the unitary matrices C/W are defined as: 
M® - i 



U!i, :lm= (^±l) . ,3!) 



I'm' dm 



This is the general relation we are looking for in the case of two-channel 
scattering. Obviously, when the off-diagonal matrix elements of the S-matrix, 
i.e. Sil and S^l vanish, Eq. (28) reduces to the famous Liischer's formula [8] 
for the single-channel elastic scattering. 



5 Eigenstates with definite cubic symmetry 



The general result (30) obtained in the previous section can be further sim- 
plified when we consider irreducible representations of the symmetry group of 
the cubic box. We know that energy eigenstates in a box can be character- 
ized by their transformation properties under the symmetry group of the box. 
For this purpose, one has to decompose the representations of the rotational 
group with angular momentum / into irreducible representations of the corre- 
sponding symmetry group of the box. For a symmetric cubic box, the relevant 
symmetry group is the cubic group O(Z). 

In a given symmetry sector, denoted by the irreducible representation T, the 
representation of the rotational group with angular momentum I is decom- 
posed into irreducible representations of O(Z). This decomposition may con- 
tain the irreducible representation T. We may pick our basis of the repre- 
sentation as: |T, a;l,n). Here a runs from 1 to dim(T), the dimension of the 
irreducible representation T. Label n runs from 1 to the total number of oc- 
currence of T in the decomposition of rotational group representation with 
angular momentum I. The matrix M is diagonal with respect to V and a by 
Schur's lemma. 

For the two-particle eigenstates in the symmetry sector T, the general for- 
mula (30) reduces to: 



= 0, 



(32) 



Here U(T) represents a linear operator in the vector space Ha(T) 4 . This 
vector space is spanned by all complex vectors whose components are V[ n , 

4 Please refer to Ref. [8] for details. 
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with I < A, and n runs from 1 to the number of occurrence of T in the 
decomposition of representation with angular momentum I [8]. To write out 
more explicit formulae, one therefore has to consider decompositions of the 
rotational group representations under appropriate cubic symmetries. 

The basic symmetry group for a cubic box is the group O(Z), which has 2 
one-dimensional (irreducible) representations A\ and A2, a two-dimensional 
irreducible representation E, and 2 three-dimensional representations T\ and 
T 2 . 5 Up to angular momentum I = 4, the representations of the rotational 
group are decomposed according to: 



= A+, 

1 = Tf, 

2 = T+ + E + , (33) 

4 = A+ + E + + T+ + T+ , 

In most lattice calculations, the symmetry sector that is easiest to investigate 
is the invariant sector: Af. We therefore focus on this particular symmetry 
sector. We see from Eq. (33) that, up to / < 4, only s-wave and g-wave 
contribute to this sector. This corresponds to two linearly independent, homo- 
geneous polynomials with degrees not more than 4 which are invariant under 
O(Z). The two basis polynomials can be identified as 3^oo and 3^40- 6 

In the first order approximation, if we neglect the mixing between the s-wave 
and g-wave, we have for the Af sector: 

det ( ^ - Voe 2 ^ ^ v / S5v /rr ^ e ^ <l)+ ' i ^ 
6 V*V^V /rr ^ M e^- Vo e™°> 




where we have also used the special parametrization (14) for the s-wave S- 
matrix elements and we have defined: 



2iAi Moo(kf)+i 



which may also be expressed as: 

cotA t = M 00 (kf) = Zoo ^ ) . (36) 

5 The notations of the irreducible representations of group O(Z) that we adopt 
here follow those in Ref. [8]. 

6 Our conventions for the spherical harmonics are taken from Ref. [22]. 
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Note that the quantities appearing above, namely 5°, S 2 , i]o, Ai and A2, are all 
functions of the energy: E = kf/(2mi) = E? + k 2 /(2m 2 ). Expanding Eq. (34) 
we get, after some algebra: 

cos(A 1 + A 2 - 5° - 5° 2 ) = Vo cos(A 1 - A 2 - 5° + 5° 2 ) . (37) 



This is the simplified formula for the s-wave ^-matrix elements. Another equiv- 
alent way of writing the same formula is: 

tan(Ai - 5?) tan(A 2 - 5° 2 ) = . (38) 



Therefore, if we neglect contaminations from higher angular momentum (mainly 
from / = 4), the parameters in the two-channel ^-matrix elements, namely 770, 
8±, 5 2 and the exact two-particle energy E satisfy a relation given by (37). 
Unlike the single channel case, where the ^-matrix has only one parameter 
(phase shift) and it is related to the exact energy in a one-to-one fashion, the 
two channel ^-matrix now has 3 real parameters and these parameters are 
related to the exact energy E by one relation. This relation is helpful since 
it provides a constraint on the four physical quantities. For example, if we 
can measure the exact energy E in lattice calculations, and if we know the 
values of 5° and S 2 from experimental data (e.g. by partial wave analysis), we 
in principle can infer information about the parameter i] which is difficult 
to measure in the experiment. If the experimental information is inadequate, 
say both S 2 and r/o are poorly determined, our result (37) still helps to setup 
a constraint between the two poorly determined physical quantities. This is 
more or less the situation in tttt scattering just above the KK threshold. Note 
that above the two-particle inelastic threshold, the number of states within 
a particular energy interval is roughly twice as many as in the single-channel 
case and every energy eigenvalue satisfies Eq. (37). 

It is instructive to discuss the situation just above the inelastic threshold. 
Assuming that the inelastic scattering only occurs in the s-wave, the total 
reaction cross section just above the threshold is given by: 

a l =° ~ A\J E - E T . (39) 

for very small (E — Ejn) > where A is some proportionality constant. This 
means that the parameter 770 behaves like: 

Vo^l-^^^E~ T , (40) 

7T 



just above the threshold. On the other hand, since the physical quantity 5 2 
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vanishes at the threshold, we may parameterize it as: 



tan So — koa 



(2) 



2«0 



(41) 



(2) 

Quantity a might be called the scattering length in the second channel. 
Substituting these into our general formula (38), we get: 7 



,(2) 



a n = 



itiiEtA 



cot 



S^Et) - A 1 (E T ) 



(42) 



If we consider the mixing of the g-wave, the formula obtained above becomes 
more complicated. We can write out the four-dimensional reduced matrix 
A4(Af) whose matrix elements are denoted as: AA(Ai)u> = mu> = mm, with / 
and V takes values in and 4, respectively. Using the general formula (26), it 
is straightforward to work out these reduced matrix elements in terms of ma- 
trix elements M.i m ;i'm' ■ We find that, in the case of cubic symmetry, Eq. (32) 
becomes: 



u {l) 

"00 



o(0) 
°11 



fcimi c(0) 
Oi 



„(1) 
"04 



'12 



U, 



k21Tl2 g(°) 
fcimi 21 

(2) o(0) 
Do 



00 



^22 



(2) 
M 04 



„(1) 
"04 





M44 s 
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fcimi c( 4 ) 
Oi 



fc^^l2 M 44 





„(2) 
"04 

fcimi 21 

(2) _ o(4) 
22 



. 



(43) 



Here matrix [/ is defined as in Eq. (31) and u lv are the corresponding matrix 
elements. At this level, so many parameters enter the relation and it seems 
that the formula is useful only when other information is available. 

We also would like to remark that, although we worked out the formulae in 
a cubic box, similar relations can also be obtained for general rectangular 
box following the strategies outlined in Ref. [20,21]. Such topologies might be 
useful for the calculations of scattering phases. It is also clear that the results 
obtained in this paper can easily generalized to more than two channels as long 
as the initial and final states in the scattering are still two-particle states. 

Finally, let us speculate about possible extension to the case of massive field 
theory. In the case of single channel scattering, it was shown in Ref. [8] that the 
result obtained in the quantum-mechanical model can be carried over literally 
to the case of massive quantum field theory as long as the non-relativistic 
dispersion relations of the particles are replaced by relativistic ones and the 
polarization effects and other effects (exponentially small) are small enough. 

7 Note that A2 ~ g^/i^oo^ Q2) ~ ^2 f° r small k%, so tan(A2 — $>) ~ — tanJ^- 
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For the case of multi-channel scattering, we expect that a similar conclusion 
to hold although a proof is still lacking. If this turned out to be true, it means 
that the results obtained in this paper can also be generalized to the case of 
massive field theory apart from corrections that are exponentially small in the 
large volume limit. 

In the case of field theory, another complication arises since in the framework 
of field theory, particle numbers are not conserved. One therefore has to spec- 
ify what one means by two-particle states [8] in a finite volume. Generally 
speaking, two-particle states for elastic scattering are the discrete (but quasi- 
continuum) spectrum states above the lowest two-particle threshold. When the 
energy is increased, we may encounter another threshold where the inelastic 
channel is opened. In the case of pion-pion scattering, the lowest two-particle 
threshold is the two-pion threshold. The next lowest threshold is the four-pion 
threshold which is below the two-kaon threshold. It is clear that the formulae 
obtained in this paper would be applicable only when the coupling between 
the four pion states and the two particle (two pion and two kaon) states be- 
comes negligible. If this were the case, then when we calculate the correlation 
function matrix among appropriate two-particle operators, two-particle states 
dominate the correlation function matrix. The energy eigenstates thus ob- 
tained are also mainly composed of two-particle states. As we have said in 
the introduction, the weakness of interaction of four pion states with two pion 
states can be seen from the chiral lagrangian in which they are coupled via 
derivative couplings. In a lattice QCD simulation, the validity of this assump- 
tion might also be checked numerically by investigating the volume depen- 
dence of the energy eigenstates obtained from the corresponding correlation 
functions. This is due to the fact that the density of state for single particle, 
two-particle and four-particle states have rather different volume dependence. 
Using this technique, for example, the authors in Ref. [23] were able to argue 
that the energy eigenstates they obtained are in fact KN scattering states 
(two-particle states) and not a single-particle penta-quark state. Similar tech- 
nique can in principle be applied to distinguish the two-particle states from 
the four-particle states in a finite volume. Of course the feasibility of this can 
only be checked in a real numerical simulation. Here we can only point out 
this possibility. 



6 Conclusions 

In this paper, we have studied two-particle two-channel scattering states in 
a cubic box with periodic boundary conditions. Assuming that energy eigen- 
states are only two-particle states, the relation of the exact energy eigenval- 
ues in the box and the physical parameters in the coupled channel ^-matrix 
elements in the continuum is found. This formula can be viewed as a general- 
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ization of the well-known Liischer's formula to the coupled channel situation 
(inelastic scattering). In particular, we show that the two-channel S'-matrix 
elements in the s-wave are related to the energy of the two-particle system by 
a simple identity, if contaminations from higher angular momentum sectors 
are neglected. This relation is non-perturbative in nature and it will help us to 
establish connections between the ^-matrix parameters in the multi-channel 
scattering with the energy eigenvalues which are in principle accessible in lat- 
tice calculations. 
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A Appendix A 



In this appendix, the proof of the theorem concerning the structure of the solu- 
tion to the radial Schrodinger equation is provided. In particular, we consider 
the coupled differential equations: 



where *&(z) is the two-component wave-function and the matrix Q(z) is given 
by: 



Q(*) = — ~ + 



1(1 + 1) I 2m 1 [E-V 1 (z)\ 2miA(z) 



z 2 



2m 2 A*{z) 2m 2 [E' - V 2 (z)} j 



(A.2) 



We would like to study the structure of the solutions to Eq. (A.l) near z — 0. 
Using standard transformation: 

*(z) = z-^iz) , (A.3) 

we find that Eq. (A.l) reduces to: 
d 2 <S>(z) ld$(z 



dz 2 z dz 



+ Q(z)$(z) = , (A.4) 
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with Q(z) given by: 



QW= J' + f) 2 + ^dE-V^)] WW 

2m 2 A*(z) 2m 2 [£' - V 2 (z)\ 



We now proceed to find canonical solutions of the type: 



= E X»* B , (A.6) 

ra=0 



where Xn are two-component coefficients, p is the index of the canonical so- 
lution. Assuming z 2 V(z) is analytic near z = and lim 2 _ > o[- 2;2 ^ / (- 2 )] = 0, we 
have: 

oo 

z 2 Q{z) = Y d QnZ n , Q = -(/ + l/2) 2 (A.7) 

n=0 



Comparing the coefficients of each order we get: 

n 

\{p + nf) - (I + 1/2) 2 ] X „ + E QfeXn-fe = . (A.8) 

fe=i 



In particular, when n = the above equation gives the index equation satisfied 
by p: 

p 2 -(l + 1/2) 2 = . (A.9) 



This equation yields two solutions: p = ±(/ + 1/2). Note that Xo is non- 
vanishing, this shows that the solutions near z = are categorized into two 
classes. In one class, the solutions behave like: <&(z) ~ z l+l l 2 ; in the other 
class the solutions behave like: <&(z) ~ z~^ l+1 ^\ Or in terms of the solution 
^(z), these two cases behave like ^(z) ~ z l and *&(z) ~ z~( l+1 ^ respectively. 
Also note that when the index p take the values of ±(Z + 1/2), the matrix 
Qo becomes identically zero. Obviously, we can find two linearly independent 
solutions Xo f° r eac h p- Therefore, we have shown that the canonical solutions 
indeed have the properties as indicated in the theorem. Since the coefficients 
Xn are fully determined by the recursion relation (A.8), all what remains to 
be shown is that the series solution (A.6) converges uniformly and absolutely 
in the neighborhood of z — 0. 
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To show the convergence of the series solution, we see from Eq. (A. 8) that: 
1 n 

Xn = =F^ TI yi;Q fc Xn- fc . (A. 10) 



Note that z 2 Q(z) is analytic near z — 0, therefore we can always find two real 
numbers M > 1 and i? > such that: 

\\Q k \\<MR- k , (A.ll) 



where || • || indicates a matrix norm. 8 Let us assume that: 

\ X u\ <M"R-»\ Xo \, i/ = l,2,.-.,(n-l). (A.12) 

we now show that the above inequality then holds for v = n. This is seen by: 



1 n 

IXn| -n(2/ + l)£ l|gfe|IIXn - fc| 

I I D— n n 

^ 7o L-n E M "~ fc+1 ^ Mn ^" n |Xo| • (A.13) 



fc=i 



Therefore, by induction, the inequality |x„| < M n i? _n |x | is true for any inte- 
ger n. It then follows trivially that the series (A. 6) is absolutely and uniformly 
convergent in a small neighborhood around z = which completes our proof 
of the theorem. 
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